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ABSTRACT 

We analyze the anisotropic clustering of massive galaxies from the Sloan Digital Sky Survey 
III Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 9 (DR9) sample, which 
consists of 264283 galaxies in the redshift range 0.43 < z < 0.7 spanning 3275 square de- 
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grees. Both peculiar velocities and errors in the assumed redshift-distance relation ("Alcock- 
Paczynski effect") generate correlations between clustering amplitude and orientation with 
respect to the line-of-sight. Together with the sharp baryon acoustic oscillation (BAO) stan- 
dard ruler, our measurements of the broadband shape of the monopole and quadrupole corre- 
lation functions simultaneously constrain the comoving angular diameter distance (2190 ± 61 
Mpc) to z = 0.57, the Hubble expansion rate at z = 0.57 (92.4 + 4.5 km s"' Mpc"'), and the 
growth rate of structure at that same redshift (dcrg/dlna = 0.43 + 0.069). Our analysis pro- 
vides the best current direct determination of both Da and H in galaxy clustering data using 
this technique. If we further assume a ACDM expansion history, our growth constraint tight- 
ens to dcTg/dlna = 0.415 + 0.034. In combination with the cosmic microwave background, 
our measurements of Da, H, and dcrg/dlna all separately require dark energy at z > 0.57, 
and when combined imply Qa = .74 ± 0.016, independ ent of the Universe's evolution at 
z < 0.57. In our companion paper (ISamushia et al.ll2012h . we explore further cosmological 
implications of these observations. 

Key words: cosmology: large-scale structure of Universe, cosmological parameters, galaxies: 
haloes, statistics 



1 INTRODUCTION 



Measurem ents of the cosmic distance-redshift rela t ion using su- 
pemo v ae JRiess et al.l 1 19981: iPerlmutter et al] Il999l : iKessler et al] 
I2OO9I : lAmanullah et al. |2010|) . the cosmic micro wave back- 
groun d JLarson et al.l uOllh . the Hubble constant (Riess e t al. 
20111) . and baryon aco u stic o scillatio ns (BAO; [Eisenstein et al 



20051: ICole et al.l l2005l: iHutsil l200d: IPadmanabhan et alj [200 



Blake etalj l2007l: iPercival et alj l2007l: lOkumura et al.l I2OO8I: 
Gaztanaga. Cabre & Huil l2009l: iKazin et alJ l20ld: IPercival et all 



2OIOI : iReid et alj I2OI0I : felake. et alj I2OI 1 J : iBeutler et alj I2OI ll) 
have revealed that the expansion of the universe is accelerating; 
either the energy density of the universe is dominated by an ex- 
otic "dark energy", or general relativity requires modification. The 
observed anisotropic clustering of galaxies can help distinguish be- 
tween these possibilities by allowing simultaneous measurements 
of both the geometry of the Universe and the growth rate of struc- 
ture. 

Galaxy redshift surveys provide a powerful measure- 
ment o f the growth rate through redshift-space distortions 
(RSD) jKaisea 119871 ^. Although we expect the clustering 
of galaxies in real space to have no preferred direction, 
galaxy maps produced by estimating distances from red- 
shifts obtained in spectroscopic surveys reveal an anisotropic 
galaxy distribution JCole, Fisher. & Weinberg||l995 : Peacock et al 



200 ll: IPercival et alj|20M: Ida Angela et alj|2008l: lOkumura et al 
2008"; 'Guzzo et al."2008l: ISamushia. Percival. & Raccanellil I20TII ; 
Blake, et al. 2011b). This anisotropy arises because the reces- 
sion velocities of galaxies, from which distances are inferred, in- 
clude components from both the Hubble flow and fr om pecu- 
liar v elocities driven by the clustering of matter (see iHamiltonl 
1998 1 for a review). Despite the fact that galaxy light does not 
faithfully trace the mass, even on large scales, galaxies are 
expected to act nearly as test particles within the cosmolog- 
ical matter flow. Thus the motions of galaxies carry an im- 
print of the rate of growth of large-scale structure and al- 
low us t o both pro be dark energy and test General Relativity 
(see e.g. I Jain & Zhang 2008; Neseris & Perivolaropoulosl 2008); 
Song & Kova md l2009l: ISong & Per cival 2009 ;: IPercival & Whitd 
2008; M cDonald & SeliakI 12 009'n'white . Song & Percival! l2009l : 
Song et all I2OI0I; IZhao et alj "|2010,: Song et aHl201lL for recent 
studies). 



The observed BAO feature in the power-spectrum and cor- 
relation function of galaxies has been used to provide strong 
constraints on the geometry of the Universe. While the BAO 
method is exp ected to be highly robust t o systematic uncertain- 
ties (see, e.^., [Eisenstein & White l2004l : IPadmanabhan & Whitd 
I2OO9I : ISeo et alj|2010l : iMehta et al.lf2011l) , it does not exploit the 
full information about the cosmological model available in the 
two-dimensional clustering of galaxies. Additional geometric in- 
formation is available by comparing the amplitude of clustering 
along and perpendicular to the line-of-sight (L OS); this is known 
as the Alcock-Paczynski (AP) test ( JAlcock & Paczvnskil 1 19791 ; 
iBallinger. Peacock. & Heavens! 1 1996h . RSD and AP tests rely on 
the measurements of anisotropy in the statistical properties of the 
galaxy distribution and are partially degenerate with each other, so 
that constraints on the growth of structure from RSD depend on 
the assumptions abou t the background geometry and vice versa 
JSamushia et alj|201 ll) . However, given high precision clustering 
measurements over a wide range of scales, this degeneracy can be 
broken since RSD and AP have different scale-dependences. Re- 
cently, the WiggleZ survey (Drinkwater et al. 2010) has performed 
a joint RSD and AP analysis that constrains the expansi on history in 
4 bins across 0.1 < z < 0.9 at the 10-15 p er cent level JBlake. et alj 
l201ld) . Using the SDSS-11 LRG sample, IChuang & Wand J2OI l!) 
perform a similar analysis to the present work. They measure the 
angular diameter distance /^^(z = 0.35) = 1048^jg Mpc and the 
Hubble expansion rate H(z = 0.35) = 82.1+^^ km s"' Mpc"' at 
z = 0.35 after marginalizing over redshift space distortions and 
other parameter uncertainties. 

Obtaining reliable cosmological constraints from the RSD and 
AP effects demands precise modeling of the nonlinear evolution 
of both the matter density and velocity fields, as well as the ways 
in which the observed galaxies trace those fields. The halo model 
of large-scale structure and variants thereof assume that galaxies 
form and evolve in the potential wells of dark matter halos, and 
provides a successful means of parametrizing the relation between 
the galaxy and halo density and velocity fields. Our modeling ap- 
proach, based on Reid & White ( 2011) , uses perturbation theory to 
account for the nonlinear redshift space clustering of halos in the 
quasilinear regime as a function of cosmological parameters, and 
then uses the halo model framework to motivate our choice of nui- 
sance parameters describing the galaxy-halo relation. We test these 
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assumptions with a large volume of mock galaxy catalogs derived 
from A'-body simulations. 

Th e ongoing Baryon Oscillation Spectroscopic Survey 
(BOSS:'Schle gel. White & Eisenstein"20 09|), which is par t of Sloan 
Digital Sky Survey III (SDSS-III) lEisenstein et alj|201lh . is mea- 
suring spectroscopic redshifts of 1.5 million galaxies, approxi- 
mately volume limited to z ^ 0.6 (in addition to spectra of 150,000 
quasars and various ancillary observations). The galaxies are se- 
lected from the multi-color SDSS imaging to probe large-scale 
structure at intermediate redshift; they trace a large cosmological 
volume while having high enough number density to ensure that 
shot- noise is not a dominant contributor to the clustering variance 
( I Whi te et al. 2011). The resulting clustering measurements provide 
strong constraints on the parameters of standard cosmological mod- 
els. 

We use the CMASS sample of BOSS galaxies that will be in- 
cluded in SDSS Data Release 9 (DR9) to constrain the growth of 
structure and geometry of the Universe. We apply RSD and AP tests 
to the data to measure the growth rate, the Hubble expansion rate, 
and the comoving angular diameter distance at Zeff = 0.57. We im- 
prove on previous, similar analyses in a number of ways. First, we 
use a model for the nonlinear anisotropic correlation function that 
is accurate to well below our statistical errors over the wide range 
of scales between 25 A^'Mpc and 160 A^'Mpc, which we validate 
using 68(/!"'Gpc)' of A'-body simulations populated with mock 
CMASS galaxies. Next, rather than assuming a fixed underlying 
linear matter power spectrum, we use a prior on P(k) based on the 
WMAP7 cosmic microwav e background constraints (ILarson et al.l 
I2OI iMKomatsu et alj201 ll) and marginalize over the remaining un- 
certainties for all fits. In addition to joint constraints on the geo- 
metric (BAO and AP) and peculiar velocity (RSD) parameters, we 
present three "null" tests of the ACDM model. To begin, we sim- 
ply ask whether any points in the ACDM parameter space allowed 
by WMAP7 provide a good fit to the CMASS clustering; in this 
case, the only free parameters are those describing how galaxies 
trace matter. In the other two tests, we fit for the amplitude of pe- 
culiar velocities using WMAP7 priors on geometric quantities, or 
we fit for the geometric parameters with WMAP7 priors on RSD. 
Thus we can present the statistical precision with which our data 
measure either peculiar velocities or the AP effect in the ACDM 
model. Finally, given the strong detection of the BAO feature in the 
monopole of the correlation function, we can break the degeneracy 
between (1 + Zeff)£>A and H with our AP measurement. We present 
the most constraining measurement of //(Zeff) from galaxy cluster- 
ing data to date using this technique, even after marginalizing over 
the amplitude of the RSD effect. As our constraints exploit the full 
shape of the monopole and quadrupole correlation functions, they 
rely on further assumptions about the cosmological model: Gaus- 
sian, adiabatic, power-law primordial perturbations, the standard 
n umber A'etr = 3.04 of massless neutrino species (see discussion 
in lKomatsu et alj|201 11 ). and that dark matter is "cold" on the rela- 
tively large scales of interest. 

This paper is organized as follows. In Sec. |2] we describe 
BOSS DR9 CMASS data and in Sec. [3] we describe the measure- 
ments of two-dimensional clustering of galaxies used in this anal- 
ysis. Sec. [previews the theoiy of the RSD and AP effects, and de- 
scribes the theoretical model used to fit our measurements. Sec.|5] 
outlines the methods we use to analyse the data and Sec.|6]presents 
the results of analysis. We conclude by discussing the cosmological 
implications of our results in Sec.|7] 




Figure 1. Number density as a function of redshift for the CMASS galax- 
ies used in this analysis. After accounting for our weighting scheme, the 
effective redshift for galaxy pairs in this sample is z^ff = 0.57. 



2 DATA 

BOSS targets for spectroscopy lu minous galaxies selected from the 
multi-color SDSS ima ging (Fuk ugita et al.ll996llGunn et al.ll99"3 : 
lYork et al.l200(ll : lGuim et al.i2006n. The targe t sele ction algorithms 
are su mmarized by lEisenstein et al.l (I2OIII) and I Anderson et al.l 
((2012). For the galaxy sample referred to as "CMASS", color- 
magnitude cuts are applied to select a ro ughly volume- l imited sam- 
ple of massive, luminous galaxies; see lMasters et al.l ( l201lh for a 
detailed examination of the properties of CMASS targets in the 
COSMOS field. The majority of the galaxies are old stellar systems 
whose prominent 4000 A break in their spectral energy distribu- 
tions makes them relatively easy to select in multi-color data. Most 
CMASS galaxies are central galaxies residing in dark matter halos 
of 1 ' ' A" ' Mq , but a non-negligible fraction are satellites which live 
primarily in halos about 10 times more massive (White et al. 201 ij). 
These galaxies are intrinsically v ery luminous and at t he high mass 
end of the stellar mass function i JMaraston et ahllprepl). Galaxies i n 



the CMASS sample are highly biased (b ~ 2. IWhite et al.ll201lh . 
In addition, they trace a large cosmological volume while having 
high enough number density to ensure that shot-noise is not a dom- 
inant contributor to the clustering variance, which makes them par- 
ticularly powerful for probing statistical properties of large-scale 

stru cture. 

lAnderson et al.l ( l2012h details the steps for generating the 
large-scale structure catalog and mask for DR9, which includes the 
data taken by BOSS through luly 2011 and covers 3275 deg^ of 
sky. In our analysis we use galaxies from the BOSS CMASS DR9 
catalog in the redshift range of 0.43 < z < 0.70. The sample in- 
cludes a total of 264,283 galaxies, 207,246 in the north and 57,037 
in the south Galactic hemispheres. Figure[T]shows the redshift dis- 
tribution and Figure|2|shows the angular distribution of the galaxies 
in our sample. 



3 MEASUREMENTS 

3.1 Two point statistics 

To compute redshift space separations for each pair of galaxies 
given their angular coordinates and redshifts, we must adopt a cos- 
mological model. We use the same one as used to generate our 
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Figure 2. BOSS DR9 sky coverage. The light gray region shows the ex- 
pected total footprint of the survey, while the colors indicate completeness 
in the DR9 spectroscopy for regions included in our analysis. Dark grey 
regions are remove d from the analysis by completeness or redshift failure 
cuts; see Sec. 3.5 of lAnderson et al.l J2012h for further details. 



mock catalogs, namely spatially-flat ACDM with fln, = 0.274 . 
This is the same cosmology as assumed in Iwhite et al. j201lh 
and in our com panion papers (Anderson et al. 2012; Manera et al. 
20121; iRossetal. 2012 ; Sanchez et al. 2012 ; Samushia et al. 20 1 2 ; 
Toieiro et al.ll2012n . Our model accounts for this assumption and 
scales the theory prediction accordingly when testing a cosmologi- 
cal model with a different distanc e-redshift relation; see Sec. 14. 51 

Using the galaxy catalog of lAnderson et al.l ( 120121) we com- 
pute weighted "data-data" (DD) pair counts in bins of s and jx 



Nm A'lm 

DD{si,fij) = 22 &k.i{Si,JXj)wt(Si)wi(Si), 

k=l l=k+l 



(1) 



where s is the comoving pair separation in redshift space, /u is the 
cosine of the angle between the pair separation vector and the LOS 
and wt is the weight of A:* galaxy in the catalog. The double sum 
runs over all galaxies and &i,i(si,ijj) = 1 if the pair separation be- 
tween two galaxies falls into bin Sj,iUj, and is zero otherwise. 

Three distinct effects contribute to the final weight vv, of 
each galax y. These weights are described in more detail in 
[Anderson et al. (2012) and Ross et al. (2012) . First, galaxies lack- 
ing a redshift due to fiber collisions or because their spectrum was 
not adequate to secure a redshift are accounted for by upweighting 
the nearest galaxy by weight w = (1 + n), where n is the number 
of near neighbors wit hout a redshift. S econd, we use the minimum 
variance weighting of lHamiltonI ( 119931) , 



w(s) = 



1 



(2) 



l+J3(s)h(z) 

where n(z) is the expected number density of galaxies at the redshift 



Ms) = 2;r 



Jo 



'ds'dii ^(s',^) 



(3) 



is the angularly averaged redshift-space correlation function inte- 
grated up to the separation of galaxies in the pair. For every galaxy 
this weight will vary depending on which pair-counting bin it is 
assigned. For a constant radial selection function this weighting 
scheme results in the minim al variance of th e estimated correla- 
tion fu nction (for details seelHamiltonlll993l). Note tha t the anal- 
yses in lSanchez et al.l ( 120 12h and lAnderson et al.l ( 120 12h use scale- 
independent weights; differences between the approaches are small 
in practice. 

The third weight corrects for angular systematics, related to 
the angular variations in density of stars that make detection of 
galaxies ha rder in areas of s ky closer to the Galactic equator (for 
details see IRoss et al.l 120121) . The total weight is the product of 
these three weights. We bin i in 23 equal logarithmic bins between 
iniin = 25.1 and Sy^.^ = 160 /i^'Mpc with dlogio.? = 0.035, and 200 
equally-spaced yu bins between and 1 . We compute "data-random" 
(DR) and "random-random" (RR) pair counts as in Eqn. \T\ except 
that each random point is assigned only the Jjis) weight and not the 
close-pair correction and systematic weights. Positions of objects in 
our random catalog are generated using the observational mask and 
redshifts are generated by picking a random redshift drawn from the 
measured redshifts of observed galaxies. Our random catalogs con- 
tain approximat ely 70 times mor e objec ts than the galaxy catalog. 

Following iLandv & Szala\i ( Il993h . the pair counts are com- 
bined to estimate the anisotropic correlation function as: 



iLsiSi,!^]) = 



DD(si,iJj) - 2DR(Si,tij) + RR(s„iij) 
RR(s„^ij) 



(4) 



Figure [3] shows our measurement of ^is(si,fij) in terms of LOS 
separation r„ = s^s and transverse separation r^ = s(l - yu^)''-. 
The central "squashing" is due to peculiar velocities. In the left 
panel, the BAO ridge at ~ 100 h^' Mpc is evident. In the right 
panel, we show the clustering signal on smaller scales; the "finger- 
of-God" effect is visible for small transverse separations but small 
on the scales we analyse. The innermost contour in the right panel 
indicates the value of fo in the smallest separation bin included in 
our cosmological analysis. 

Rather than work with the two-dimensional correlation func- 
tion ^(.s,yu,), we conduct our cosmological analysis on the first two 
even Legendre polynomial moments, ^o(^) and <f2(.5), defined by 

2f+ 1 



&("■) 



I 



djJ, ^(i',yUj)Lf(//.,), 



or equivalently. 



^{s,ti,) = Y,^e(s)Lf(M,). 



(5) 



(6) 



Here Lf is the Legendre polynomial of order L By symmetry all 
odd-^ moments vanish and on large scales the measurements be- 
come increasingly noisy to larger £. The correlation functions iois) 
and i2(s) are estimated from ^isi-^hl^j) using a Riemann sum to 
approximate Eqn. |5] We include all galaxy pairs between 25 and 
160 /z^'Mpc in our analysis. We also caution the reader that we have 
adopted logarithmical l y-spaced bins, while our companion p apers 
JAnderson erZll2012l: Iross et alJlIoT^ ; ISanchez et alj|2012h ana- 
lyze clustering in linearly- spaced bins of differing bin sizes. Our 
measurements of ^o and ^2, along with d iagonal errors estimated 
from mock catalogs ( IManera et al.ll2012l ; see Sec. 13.2) are shown 
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Figure 3. Left panel: Two-dimensional correlation function of CMASS galaxies (color) compared with the best fit model described in Section |6?Tl (black lines). 
Contours of equal f are shown at [0.6, 0.2, 0. 1, 0.05, 0.02, 0]. Right panel: Smaller-scale two-dimensional clustering. We show model contours at [0. 14, 0.05, 
0.01, 0]. The value of fo at the minimum separation bin in our analysis is shown as the innermost contour The p x \ "finger-of-god" effects are small on the 
scales we use in this analysis. 



in Figure |4] The effective redshift of weighted pairs of galaxies in 
our sample is z = 0.57, with negligible scale dependence for the 
range of interest in this paper. For the purposes of constraining cos- 
mological models, we will interpret our measurements as being at 
z = 0.57. 



3.2 Covariance Matrices 

The matrix describing the expected covariance of our measure- 
ments of f f (.s) in bins of redshift space separation depends in linear 
theory only on the underlying linear matter power spectrum, the 
bias of the galaxies, the shot-noise (often assumed Poisson) and the 
geometry of the survey. We use 600 mock galaxy catalogs, based 
on Lagrangia n perturbation theory (LPT) and described in detail in 
anera et alj ( 12012) . to estimate the covariance matrix of our mea- 



matrix, or "reduced covariance matrix", defined as 



Im 

surements. We compute ^f(.5,) for each mock in exactly the same 
way as from the data (Sec. l3.lt and estimate the covariance matrix 
as 



^'7^' = ^ Z (^A ^■'') - ^^. (■^')) (4(-'.) - ^^.(-^v)) ' 



(7) 



where iff (.?,) is the monopole {E = 0) or quadrupole (( = 2) correla- 
tion function for pairs in the ;''' separation bin in the k''' mock. if(s) 
is the mean value over all 600 mocks. The shape and amplitude of 
the average two-dimensional correlation function computed from 
the mocks is a good match t o the measured corr elati on function 
of the CMASS galaxies; see iManera et akl ( 1201 2|) and IRoss et al.l 
(120120 for more detailed comparisons. The square roots of the di- 
agonal elements of our covariance matrix are shown as the error- 
bars accompanying our measurements in Fig. |4] We will examine 
the off-diagonal terms in the covariance matrix via the correlation 






iW' 



C"; 



(8) 



where the division sign denotes a term by term division. 

In Figure |5] we compare selected slices of our mock covari- 
ance matrix (points) to a simplified prediction from linear theory 
(solid lines) that assimies a constant number density n = 3 x lO"'* 
(/i~' Mpc)~^ an d neglects the effe cts of survey geometry (see, e.g., 
iTegmarkI Il997h . IXu et al.l (120121) performed a detailed compari- 
son of linear theory predictions with m easurements from th e Las 
Damas SDSS-II LRG mock catalogs (McBride et alj Iprecl) . and 
showed that a modified version of the linear theory covariance with 
a few extra parameters provides a good description of the W-body 
based covariances for ^q{s). The same seems to be true here as 
well. The mock catalogs show a deviation from the naive linear 
theory prediction for £,2{s) on small scales; a direct consequence is 
that our errors on quantities dependent on the quadrupole are larger 
than a simple Fisher analysis would indicate. We verify that the 
same qualitative behavior is seen for the diagonal elements of the 
quadrupole covariance matrix in our smaller set of A'-body simu- 
lations used to calibrate the model correlation function. This com- 
parison suggests that the LPT-based mocks are not underestimating 
the errors on ^2, though more A'-body simulations (and an account- 
ing of survey geometry) would be required for a detailed check of 
the LPT-based mocks. 

The lower panels of Figure |5] compare the reduced covari- 
ance matrix to linear theory, where we have scaled the O^^ pre- 
diction from linear theory down by a constant, c,. This compar- 
ison demonstrates that the scale dependences of the off-diagonal 
terms in the covariance matrix are described well by linear the- 
ory, but that the nonlinear evolution captured by the LPT mocks 
can be parametrized simply as an additional diagonal term. Finally, 
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Figure 5. Upper panels: Diagonal elements of the monopole and 
quadrapole components of the covariance matrix computed from the LPT- 
based mocks (points) compared with the linear theory prediction (soHd 
lines). Lower panels: Two slices through the reduced covaiiance matrices 
C™"'' and Cf •'■'='' for separation bins of 33 and 103 ft"' Mpc. Linear theory 
predictions for the reduced covariance matrices in the lower panels have 
been scaled by a constant factor to produce good agreement between linear 
theory and mock covariances for off-diagonal elements, demonstrating that 
the scale dependence of the off-diagonal terms matches the mock covari- 
ance matrix well, but that there is extra diagonal covariance in the mocks 
compared with linear theory. Elements of C are small (not shown). 



4 THEORY 

4.1 Redshift Space Distortions: Linear Theory 

The effects of redshift space distortions in the linear regime are 
well-known (Kaiser 1987; Fisher 199^ see also|Hamilton 1998 for 
a comprehensive review). We briefly summarize them here. The 
redshift-space position, s, of a galaxy differs from its real-space 
position, X, due to its peculiar velocity. 



Figure 4. ^o(i) and f2(-«) measured from BOSS CMASS galaxies. The er- 
rorbars correspond to diagonal elements of the covariance matrix. The best- 
fitting model described in Section 16^1 is shown as the solid curve. We use 
23 logarithmically spaced bins and include pairs between 25 and 160 ft"' 
Mpc. 



while not shown here, the reduced covariances between ^o and ^2 
are small. 

Our analysis uses the LPT mock-based covariance matrix, 
which accurately accounts for both complexities of the survey ge- 
ometry as well as nonlinear corrections to the growth of structure 
on the relatively large scales of interest here, and this allows us to 
accurately report uncertainties associated with both our measure- 
ments and parameter fits. 



s = X + v.(x)z, 



(9) 



where v-(x) = uJ(%)l(aH) is the change in the apparent LOS posi- 
tion of a galaxy due to the contribution of the LOS peculiar veloc- 
ity ;(; to the galaxy's redshift. Since overdensities on large, linear 
scales grow in a converging velocity field (V • v = -f6„,), the effect 
of peculiar velocities induces a coherent distortion in the measured 
clustering of galaxies that allows us to measure the amplitude of 
the peculiar velocity field. In linear theory, and with some approx- 
imati ons, the anisotropic galaxy power spectrum becomes (IKaiseil 
[l987h 



PliKfik) ^ 



: [b + fllff P'-Jk) = b'{\ +lillfj' P'-Jk) 



(10) 



where b is the linear galaxy bias, Sg = bS,,,, f = dlncrg/dlna is the 
logarithmic growth rate of matter fluctuations, and yjj. is the cosine 
of the angle between k and the LOS. 
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4.2 Legendre Moments of ^(r) 

In linear theory (Eqn. 1101 . only the ^ = 0, 2 and 4 moments con- 
tribute to the power spectrum P'L(k), and its Fourier transform ^'l{s). 
The two are simply related by £''' order Bessel functions: 

-^Pi{k)ic{ks). (11) 

Given a tight constraint on the underlying shape of the lin- 
ear matter power spectrum, the two-di mensional cluster i ng of 
galaxies constrains both bag and fa^, flPercival & White! l2008l : 
[White, Song & Perc ivall l2009h . In this work we measure and 
model only the monopole and quadrupole moments of the cor- 
relation function, ^o,2(-^)- These two moments are sufficient to 
constrain both berg and /erg and encompass most of the avail- 
able information on the peculi ar velocity field for t he highly bi- 
ased galaxies of interest here ( IReid & White! 1201 II) . in addition 
to b eing more easily modelled than h igher-^ moments (but see 
also !Taruva. Saito. & Nishimichi! !201l!) . Therefore, we collapse 
the anisotropic clustering information in ^{Va-, r„) into two one- 
dimensional functions ^o.2(^) before extracting cosmological infor- 
mation from the anisotropic galaxy clustering. 

4.3 Redshift space halo clustering in the quasilinear regime 

The Kaiser formula describing the linear effect of redshift space 
distortions breaks down even on quite large scales. An accurate 
model of the two-dimensional clustering of galaxies must account 
for nonlinear evolution in the real space matter density and velocity 
fields, nonlinear galaxy bias, and the nonlinear mapping between 
real and redshift space separations of pairs of galaxies. The sim- 
plest picture of galaxy formation asserts that galaxies occupy dark 
matter halo s, and so as a step tow ards understanding the clustering 
of galaxies, !Reid & White! ( 1201 ih showed that a streaming model 
where the pairwise velocity probability distribution function is ap- 
proximated as Gaussian can be used to relate real space clustering 
and pairwise velocity statistics of halos to their clustering in red- 
shift space. We will demonstrate in Sec. 14.41 that the same model 
describes the clustering of galaxies: 

dy 



l+C{r.,r,): 



/I 



l+C{r)U 



-{rn-y-lir\l{r)]- lltr\^_(r.fi) _ 



(12) 
where r^r and r,, are the redshift space transverse and LOS distances 
between two objects with respect to the observer, y is the real space 
LOS pair separation, /j = y/r, ^ is the real space galaxy correlation 
function, v^iir) is the average infall velocity of galaxies separated 
by real-space distance r, and a^^ir,^) is the rms dispersion of the 
pairwise velocity between two galaxies separated with transverse 
(LOS) real space separation r^ (y). 

fg('')5 Vi2(r) and a\r^{r,ii) are computed in the framework of 
Lagrangian (fO and standard perturbation theories (vi2, cr^2)- Only 
two nuisance parameters are necessary to describe the clustering of 
a sample of halos or galaxies in this model: b\L = b - \, the first- 
order Lagrangian host halo bias in real space, and cTp^Q, an additive, 
isotropic velocity dispersion accounting for small-scale motions of 
halos and galaxies which will be described below. Further details 
of the model, its numerical implementation, and its accuracy can be 
found in lkeid & White ( ,2011,) and AppendixJB] 



4.4 From halos to galaxies 

iReid & WhiT3 J201ll) examined the validity of Eqn. [12] only for 



halo clustering rather than galaxies, thus the model must be ex- 
tended and checked with a realistic sampl e of mock galaxies. We 
use the machinery of the halo model (see ICoorav & Shethll2002l 
for a review) to describe the galaxy density field in terms of the 
density field of the host halos. Of particular importance for model- 
ing redshift space distortions is the distinction between "central" 
and "sateUite" galaxies (jKa uffmann. White. & Guiderdonil 1 1993 : 
ICole et al.lll994l : lKravtsov et alj2004l) . When modeling an approx- 
imately mass-limited galaxy sample, the first galaxy assigned to a 
host halo is considered central, and its position and velocity are 
that of the host halo center. Satellite galaxies orbit in the potential 
well of the host halos, and so are offset in both position and ve- 
locity from the halo center. In our particular implementation of the 
halo model, satellite galaxies are randomly drawn from the dark 
matter particle members of the host halo in our simulation. As the 
virial velocities of massive halos can be large (amounting to red- 
shift space LOS separations of tens oihT^ Mpc), intrahalo velocities 
can distort the redshift space correlation function. In the limit that 
these virial motions are uncorrelated with the quasilinear velocity 
field of interest, they can be accounted for by additional convolu- 
tion along the line of sight. 

To assess the impact of virial motions on the obs erved galaxy 
cluster ing we use the mock catalogs described in IWhite et al.l 
( I2OIII) , which closely match the small-scale clustering of CMASS 
galaxies. We compute ^o,2{s) from these mocks (the average is 
shown as the error bars in Figure [6] and recompute fo,2 after ar- 
tificially setting the intrahalo velocities to (the dashed curves, 
with the one for ^0 covered by the solid line). Intrahalo velocities 
suppress the amplitude of ^2 on the smallest scales we attempt to 
model, reaching a 10 per cent correction at 25 /i^'Mpc. The reason 
for this suppression is that, on small scales, d^' /dy < and non- 
negligible. This causes a net transport of pairs to larger separations 
in redshift space. Note this is opposite to the effect of quasilinear 
peculiar velocities, which make the separation of a pair in redshift 
space on average smaller than in real space. 

We include the effect of intrahalo velocities in our model by 
including an extra convolution with a Gaussian of dispersion ctfog- 
The solid curves in Fig.|6]are then the predictions for <fo,2 with the 
best fit (Tp^Q = (3.2/z"'Mpc)^ or 21 Mpc". The model successfully 
describes the effect of intrahalo velocities on the monopole and 
quadrupole correlation functions. 

In Appendix lB4l we quantify the impact of our uncertainties in 
the halo occupation distribution (HOD) of CMASS galaxies and the 
possible breakdown of our assumption that the first galaxy assigned 
to each halo is "central" (i.e., has no intrahalo velocity dispersion) 
on the value of our nuisance parameter cTp^j,. As a result of these 
considerations, we place a uniform prior on (Tp^^ between and 40 
(Mpc)^ 



4.5 Alcock-Paczynski Effect 

Galaxy redshift surveys collect two angular coordinates and a red- 
shift for each galaxy in the sample. A fiducial cosmological model 
must be adopted to generate maps and measure clustering as a func- 
tion of comoving separations. This mapping depends on both the 
angular diameter distance and the inverse of the Hubble parameter 
at the redshift of each galaxy pair. To a good approximation the in- 
ferred galaxy clustering in a different cosmological model can be 
obtained from the fiducial one by a single rescaling of the trans- 
verse and parallel separations jPercival et al.ll2010l) . Rather than 
modify our observed galaxy clustering, we will account for the AP 
effect when we test different cosmological models by introducing 
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Figure 6. Error bars enclose the mean s^o.li^) measured from the 

IWhite et al.1 1201 11) mock galaxy catalogs. The solid line is our model fit, 

where (ri, „ has been varied to minimize the difference. The dashed line 
FoG 

shows s^2 for the mock galaxies when their intrahalo velocities are artifi- 
cially set to 0, and indicates that intrahalo velocities suppress f 2 by ai 10 
per cent on the smallest scales we are fitting. The nuisance parameter CTp^j, 
adequately describes the effect of intrahalo velocities. Dotted lines show the 
predicted ^2 when varying the Alcock-Paczynski paraineter F(z) by ± 10 per 
cent and holding Dy (and thus fo) fixed. 



two scale parameters, tt^ and a^, into the theoretical correlation 
functions we are fitting to: 



ttj. = 









(13) 
(14) 



where Da denotes the physical angular diameter distance. Here ^"™ 
is the expected two-dimensional correlation function (/ the mea- 
sured galaxy correlation function were computed assuming the true 
redshift-distance relation in the cosmology being tested, .f"'' is the 
prediction for the measured correlation function, given that galaxy 
separations were computed using the fiducial cosmology model. 
That is, cfx and a^ scale the 'true' separations to the ones calcu- 
lated using the fiducial cosmology. 

The spherically averaged correlation function, fo> is sensitive 
to the parameter combination 



Dv(z) = m+zfDi(z) 



cz \ 



1/3 



(15) 



The quadrupole of the measured cor relation function allows a 
measurement of a second combination ( Alcock & Paczvnskill 19791 : 
IPadmanabhan & Whitell2008l : lBlake. et"ani201ld) 



F(z) = (1+z)Da(z)H(z)/c, 



(16) 



thus breaking the degeneracy between (1 + Zss)Da and H. To a 
good approximation, changing Dy simply rescales the value of s in 
the predicted correlation function, while F(z) primarily affects the 
quadrupole. Figure|6]shows the effect of varying F by ± 10 per cent 
on ^2 at fixed Dy with the dotted curves. The scale-dependence of 
A^2 due to the AP effect will allow us to separate the effects of 
peculiar velocities and the AP effect. 



5 ANALYSIS 

5.1 Cosmological Model Space 

Given an underlying linear matter power spectrum shape 
Piin(fc, Zeff ), we consider the predicted galaxy clustering as a 
function of five parameters: pg.^\ = ji'cr8,/crg,crFoG,(l + 
z,:ff)DA{Zes), H{zes)]- Since the normalization of P[m(k,Zef{) (de- 
noted throughout as erg) determines the amplitude of the second- 
order perturbation theory corrections in our model, in principle we 
should be able to separately determine fecrg, /erg, and erg. In prac- 
tice, the dependence is sufficiently small and degenerate with the 
nuisance parameters bag and o"pj,Q that the degeneracy cannot be 
broken; see Appendix IB 3 1 for details. 

We have also assumed that any error in the fiducial cosmolog- 
ical model used to compute ^o,2(-5) from the CMASS galaxy catalog 
can be absorbed in a single scaling of distances, inteipreted at the 
effective redshift of the survey. If the assumed redshift dependence 
of (1 + Zeti)DA{z) and H(z) is grossly incorrect, we would expect a 
difference in the correlation functions split on redshift; we see no 
evidence for this in our tests JRoss et al.l2012h . 



5.1.1 Prior on the linear matter power spectrum from WMAP7 

While the scale dependence of galaxy clustering itself can constrain 
the shape of the linear matter power spectrum, at present (and cer- 
tainly with the imminent public release of Planck data) the con- 
straints enabled by CMB measurements are stronger. The strong 
CMB constraints mean we can use the entire linear matter power 
spectrum as a standard ruler determined by observations of the 
CMB, rather than only the BAO feature. This approach relies on 
further cosmological model assumptions that are consistent with 
the current data, but from which moderate deviations are still al- 
lowed. 

The temperature of the CMB has been measured exquisitely 
well JMather et al.lll994h . and determines the physical energy den- 
sity in radiation, n^h^. I n the minimal cosm ological model allowed 
by current observations dLarson et al.bOllb . namely a flat ACDM 
cosmology with nearly scale-invariant scalar, adiabatic, Gaussian 
fluctuations along with the three standard, nearly massless neutrino 
species, only three additional parameters determine the shape of 
the underlying linear matter power spectrum, Pii„(lc). Relative peak 
heights in the CMB determine the physical energy densities in cold 
and baryonic matter, Clc.bh^, and the overall scale-dependence of 
the CMB power spectrum determines the spectral index n, of the 
nearly scale-invariant scalar primordial fluctuations. Constraints on 
these parameters do not rely on the distance to the cosmic mi- 
crowave background, and thus are immune to the behavior of dark 
energy at lower redshifts than the last scattering surface. Moreover, 
small -scale CMB experiments JKeisler et al.l 1 20 111 ; iHlozek et al.l 
120111) now probe fluctuations on the same scales as galaxy clus- 
tering measurements, and find no compelling evidence for, e.g., a 
running of the spectral index. Allowing for running of the spec- 
tral index would degra de the CMB constr aints on the linear matter 
power spectrum (e.g., iMehta et al.ll2012h but we will not include 
this additional parameter-obviously, our constraints should be in- 
teipreted in the context of our model assumptions. 

One important extension of the minimal cosmological model 
is allowing neutrinos to have m ass; neutrino oscillation experiments 
suggest that 2 my > 0.05 eV JAbazaiian et al.ll201lh . As the uni- 
verse expands and cools, massive neutrinos become non-relativistic 
and modify the linear matter power spectrum inferred from the 
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CMB, as well as alter the expansion history as compared with the 
massless neutrino case. In this work we ignore this additional un- 
certainty in the shape of Pii„(k). Current upper bounds th at combine 
sever a l cosmological probe s find 2 m,, < 0.3 eV (e.g.. lReid et al.l 
l20ld : Ide Putter et~ai]|2012h . whic h is safely below th e detectable 
level in the DR9 CMASS sample ^ Sanchez et alj2012l) . 

As we explore the cosmological constraints available from 
our data set, we will consider a set of power-spectrum 
shapes parametrized by p^ = [ilbh- , Cl^h^ , n^] . We will 
marginalize over parameters 0^ either by importance resamplin 
JGilks. Richardson & Spiegelhalteil Ima iLewis & Bridi3 |200: 



the public Monte Carlo Markov Chains for a spatijilly flat ACDM 
model provided by the WMAP collaboration, or by approximating 
the constraints on p^ with a three-dimensional Gaussian likelihood. 
Note that the three-dimensional constraints on linear matter power 
spectrum parameters are very nearly independent of cosmological 
model extensions that change the expansion rate only for z <^ 1000 
(i.e., introducing w or fij). 

5.2 Models 

In the present paper, we interpret the anisotropic clustering of 
CMASS galaxies only in the context of the ACDM cosmology. 
By relaxing assumptions about the redshift-distance relation and/or 
the growth of structure in the ACDM model, we report the statisti- 
cal precision with which CMASS measurements constrain both the 
peculiar velocity field /crgCZcff) and geometric quantities OACZcff) 
and H(zcfi), without adopting a particular cosmological model that 
specifies how these quantities are related. We consider the follow- 
ing four models: 

• Model 1: WMAP7+CMASS flat ACDM. In the flat ACDM 
cosmology, /crg(zeff), D^dcff) and H(zcff) are all determined once 
n„, CTg, and Ho are specified. This provides a null test of our as- 
sumptions relating CMB fluctuations to predicted galaxy fluctua- 
tions. 

• Model 2 WMAP7-I-CMASS ACDM geometry, free growth. 
The analysis for this model is the same as Model 1, except that we 
consider fcrg a free parameter in the CMASS galaxy clustering fits. 

• Model 3 WMAP7-I-CMASS ACDM growth, free geometry 
The analysis for this model is the same as Model 1, except that 
we consider -DA(Zeff) and Hizcs) as free parameters in the CMASS 
galaxy clustering fits. 

• Model 4 WMAP7-I-CMASS, free growth, free geometry: In 
this model, /crgfeif), Z?A(Zcff) and //(Zcff) are all free parameters in 
the fit to the galaxy clustering data. This allows us to determine 
how well our data break the degeneracy between the RSD and AP 
effects, and to present the constraints originating from the ampli- 
tude and scale-dependence of the galaxy quadrupole ^i(s) in the 
most model independent way possible. This multivariate distribu- 
tion can be used to constrain any model that does not alter the shape 
of the linear matter power spectrum at Zas from that inferred from 
the CMB on the relatively large scales of interest here. 

In all four scenarios we allow the value of cTp^^ to vary between 
Mpc^ and 40 Mpc^ with a flat prior and marginalize over both 
CTp^Q and bcTji when deriving final results. We detail our methods for 
sampling the multi-dimensional probability distribution functions 
of interest in Models 1 through 4 in Appendix |C] 

This approach to parameter fitting allows our estimates of the 
growth of structure and geometry to be independent of many model 
assumptions and can be used to put constraints on more general 
models of gravity and dark energy. However, they still rely on the 



standard model in three ways. First, we assume that the processes 
in the early Universe that were responsible for setting up the lin- 
ear matter power spectrum at recombination do not change signifi- 
cantly, which is true for the most popular models of modified grav- 
ity and dark energy. Second, we assume that growth at the level of 
linear perturbation theory is scale independent between the CMB 
epoch and the effective redshift of our sample. Third, we use GR to 
compute the perturbation theory corrections to the galaxy cluster- 
ing predictions. The perturbation theory corrections are not large, 
and they are most important on small scales where cTp^^ also be- 
comes important; therefore we cannot strongly constrain the am- 
plitude of the higher-order corrections. Constraints on models with 
scale-dependent growth should be derived directly from the corre- 
lation function measurements and their covariance. 



5.3 The meaning of erg 

We follow the standard convention of denoting the amplitude of the 
matter power spectrum by ai, even though we restrict our analysis 
to scales s > 25 A"' Mpc, so a different weighted integral over 
P{k) that is concentrated on larger scales would more accurately 
reflect our constraints on the growth rate of matter fluctuati ons. In 
partic ular, since P{k) is well-determined for k in Mpc"' (e.g. I White! 
[200g), and the BAO scale provides a standard ruler with percent- 
level precision, our data constrain the amplitude of matter fluctua- 
tions on scales > 36 Mpc. In practice, the tight constraints on the 
shape of P{k) means that differences arising from how one specifies 
its amplitude are small when computing parameter constraints, as 
long as h is well-determined in the model. In Model 2, we sample 
power spectra from the WMAP ACDM chain, and take the tradi- 
tional value of CTg to relate the model parameter / and the reported 
constraint /erg. In Model 4, we do not specify a value of h with 
each sampled power spectrum, so we normalize the power spec- 
tra by fixing ctr, where R = ^jQl = 1 1.4 Mpc. For power spectra 
drawn from WMAP7 ACDM chains, cr^i/crg = 0.99±0.024; the off- 
set and variance between these parameters is negligible compared 
to our measurement errors on /erg. 



6 RESULTS 

In this section we present the results of fitting our analytic model 
for ^o.2(-5) to the observed galaxy correlation functions. Figure |7] 
summarizes our constraints from the shape of the observed angle- 
averaged correlation function ^o(s), while Figures [8]and|9]highlight 
our parameter constraints from the observed anisotropic galaxy 
clustering. Constraints on both the peculiar velocity amplitude and 
geometric quantities (I +Zcff)-DA(Zcff) and -tf(Zeff) are summarized in 
Table[T]for Models 2-4. 



6.1 Goodness of fit and ACDM results 

We include 23 separation bins for both £ = and f = 2 in our x~ 
analyses. In this section, we consider models with increasing num- 
bers of free parameters, and ask whether changes in x^ across the 
models indicate a preference for parameter values outside the pre- 
dicted values from WMAP7 in the ACDM model. We first fix the 
underlying power spectrum to the one assumed for all of the mock 
catalogs, and also fix o"p^Q = 21 Mpc^, the best fit value to our A'- 
body based mock galaxy catalogs. We vary the galaxy bias, and 
find a minimum ;t'^ = 45.7 for 45 degrees of freedom, demonstrat- 
ing that the mock galaxy catalogs used to validate our model and 
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Model 



Ao-g 



/o"8 



Dv [Mpc] F 



(1 + Zss)Da [Mpc] // [km s"' Mpc" 



2 
3 
4 
WMAP7 ACDM 



1.228+ 



a+0.033 
-0.032 
J 246+0.043 

1 TTO+Kotv 
-^ -0.050 



415+0034 



A 4^7+0.069 
"■^^ ' -0.063 



2076!« 
2070+« 



0.451 ±0.025 2009 ±42 



0.675^1 
0.6635+*^ 



2204: 

2190: 



44 
61 



2113+'3 



92.9!3| 
92.4+4:5 

94.2+tl 



Table 1. The median and 68.3 per cent confidence level intervals on parameters h(Ti , fa-g , absolute distance scale Dy (Eqn. llSK Alcock-Paczynski parameter F 
(Eqn. ll6K as well as derived parameters, comoving angulai' diameter distance ((1 +Zcff )Oa) and expansion rate (H). To obtain these constraints, we marginalize 
over tTp^Q and power spectrum shape parameters /?s = {il\,h-,il^h-,ns] for Models 2-4, as described in Section \52] We interpret our measurements at the 
effective redshift of our galaxy sample, ZeS = 0.57. 



compute our covariance matrix are consistent with the observed 
galaxy clustering. If we allow CTp^Q to vary as well, x^ = 42. 1 at 
""fog ~ ^^ Mpc^ . The difference indicates that we cannot expect a 
strong constraint on CTp^^ within our prior when other cosmologi- 
cal parameters are varying; it is important, however, to marginalize 
over this nuisance parameter, since it increases our uncertainty in 
fo-g; see the discussion in Section [631 

If we restrict ourselves to ACDIVI models consistent with 
WMAP7 (Model 1), we find a minimum ^^ value of 39.3 at n„h^ = 
0.1395 and /fo = 68.0kms"']VIpc"' when fccrg and cTp^Q are varied. 
IVIodel 4 has the most free parameters: five describing the galaxy 
clustering and three specifying the linear matter power spectrum. 
In this case we find a minimum ;^'^ of 39.0 for 41 degrees of free- 
dom. This best-fitting model is shown with our measurements of 
the correlation function in Figures |3]and|4] and has parameter val- 
ues bcTg = 1.235, /o-g = 0.437, o-p„G = 40 Mpc-, Da = 2184 
Mpc, H = 91.5kms-iMpc-', n„,h- = 0.1364, n^h- = 0.02271, 
n, = 0.967. We conclude that the observed ^0,2 is fully consis- 
tent with the ACDM cosmology; changes in;^'^ values between the 
models do not indicate a significant preference for parameter values 
of /(Tg, F, and Dy/r, outside of the values predicted by WMAP7 
in the ACDM model. 

CM ASS measurements also improve constraints compared to 
WMAP7 in the ACDM model: n„,h- = 0.1363 ± 0.0035, fl,,, = 
0.283 ± 0.017, and Ho = (69.3 ± 1.5)kms-'Mpc"'; WMAP7 alone 
finds n,„h- = 0.1334 ± 0.0056, fi„, = 0.266 ± 0.029, and Ho = 
(7 1.0 ± 2.5)kms~'Mpc~ '. Comparison with the BAO-only results 
of I Anderson et al.l 1 120121) demonstrates that in this minimal model, 
nearly all of the additional information on these three parameters is 
coming from the BAO feature. 

However, the shape of the measured galaxy correlation func- 
tion does provide an independent probe of the underlying linear 
matter power spectrum. With a strong prior on Q.i,h~ and n, taken 
from the CMB, the clustering of galaxies is sensitive to the peak 
in the linear matter power spectrum, which depends on the horizon 
size at matter-radiation equality, oc Q.i„h^ at fixed effective num- 
ber of relativistic species (A'eff = 3.04 for the standard three neu- 
trino species). The scale at which the peak appears depends on the 
low redshift distance relation, so the broadband shape of the angle- 
averaged galaxy power spectrum or correlation function constrains 
Q.„h-Dv(Zsff), in addition to the constraint on £'v(Zeff)/''s(Zdras) that 
comes from the location of the BAO feature. We illustrate the 
constraining power of our dataset by fixing ili^h^ = 0.02258 and 
n, = 0.963, and computing the CMASS-only likelihood in the 
Dv - ii„,h^ plane shown in Figure|7]as the solid contours. For this 
exercise we use only the monopole {£ = 0) measurements, and find 
a minimum ;^'^ value of 18.2 for 19 degrees of freedom (DOF). 

For comparison, we also isolate the broadb and shape infor- 
matio n by fitting to a no-wiggle power spectrum (lEisenstein & Hij 
Il998l) . which should primarily be sensitive to il,„h^Dv(,Zss)- Re- 




n 



Figure 7. Contours of A^^ = 2.30 and 6.17 for fixed fij/i^ = 0.02258 
and «s = 0.963 for the monopole {£ = 0) galaxy clustering measurements 
alone (solid). For compa rison, we also comp ute x^ using the "no-wiggles" 
power spectrum from Eisenstein & Hull ! 19981) (dash-dot) to isolate informa- 
tion from the broadband shape of the correlation function without the BAO 
feature; this fit is primarily sensitive to the apparent location of the peak 
in P(k), which corresponds to the horizon size at matter-radiation equality, 
oc n„,/i^Di/(Zcfr). We also project the WMAP ACDM constraints onto these 
parameters, and show 68 and 95 per cent contours (dashed). 



suits of this fit are shown as dash-dot contours in Figure |7] This 
model provides a poor fit to the measured correlation function, with 
xl,j„ = 38.5 for 19 DOF, indicating a strong preference for models 
with the expected BAO feature. However, we do find that the in- 
ferred value of n,„h- from the broadband shape of the measured 
galaxy correlation function is consistent with the prediction from 
the CMB. 

Projecting the with-BAO model fits onto Q.„,h^, we find 
Q.„lf = 0.142 ± 0.01 1. Translating the results of a similar analysis 
from Reid et al.l (2010) for the SDSS-II LRG sample (z^ff = 0.31) 
to the same assumptions yields il^h^ = O.Hl+Q-jjjj, while the 
WiggleZ analysis of emission line galax ies at ZeS ~ 0.6 finds 
n,„h^ = 0.127±0.011 JBlake. etal.l2011ah . Strictly speaking, these 
constraints are not uncorrelated since they have a small amount of 
overlapping volume; neglecting their small correlation, the com- 
bined galaxy clustering estimate for Q^/z^ is 0.137 ± 0.0064, which 
marginalizes over the low redshift distance-redshift relation and 
is in excellent agreement with the WMAP7 ACDM constraint of 
nji^ = 0.1334 ± 0.0056 (dashed contours in Figure|7). 
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Figure 8. One-dimensional constraints on f<T%, (1 + Zcff)OA(Zeff)> and //(Zcff) under different model assumptions. The dashed curves indicate WMAP7- 
only ACDM. The solid (models 2 and 3) and dotted (model 4) curves are constraints derived from the CMASS fo,2 measurements with a WMAP7 prior 
on the underlying linear matter power spectrum P(^/Mpc"'). The solid curves additionally use the ACDM parameters in the WMAP7 chains to fix either 
(1 + Zi;ff)OA(Zcff) and //(zeff) (left panel), or /erg (right two panels). All three constraints degrade when fitting for geometry and growth simultaneously using 
the CMASS observations. 



6.2 Constraints on the peculiar velocity field amplitude 

In Figures [8] and |9] we compare our constraints on the peculiar ve- 
locity field amplitude, /crgCzeff). from Models 2 (solid) and 4 (dot- 
ted) to the predicted distribution from WMAP7 (dashed), assuming 
a flat ACDM cosmology. The 68 per cent confidence intervals for 
Models 2 and 4 are listed in Table [T] These measurements agree 
with the ACDM WMAP7 expectation, 0.451 ± 0.025. 



6.3 Geometric constraints 

Our tightest geometric constraint comes from the BAO fea- 
ture in the monopole correlation function. The cosmologi- 
cal parameter dependence of the location of the BAO fea- 
ture is given by the sound horizon at the dr a g epo ch, rj(zdrag) 
(we use the definition in Eisenstein & HJ 119981) . We find 
(Dy(ZefF)/r,(z<irag)/(£'v(z.ff)/'-.(Zdmg))fiduciai = 1.023±0.019. The dif- 
ference between our best fit value and the p re-reconstruction fits to 
mo nopole correlat i on fun ction presented in [Anderson et al.l ( 120 12n 
and lSanchez et alj 1 120121) is due to our different choice of binning 
rather than our fitting methodology; we verified that with the same 
measurement values and covariance matrix, our method recovers 
the sa me central value as the result presented in [Anderson et alJ 
( I2OI2I) . 

Though our fits include information from the broadband 
shape of the correlation function, the resulting central value and 
error on Dylrs i zird^) are consistent with the fits performed in 
[Anderson et alJ ( |2012|) . which marginalize over the broadband 
shape of correlation function or power spectrum. We therefore con- 
clude that essentially all of the information on the distance scale 
Dv is coming from the BAO feature in the correlation function, as 
was also true in the analysis of the SDSS-Il LRG power spectrum 
JReid et alJl201Q) . This can be seen in Fig.[7] where at fixed Q.,„h^, 
the constraint on Dy is 2.5 times weaker for the "no-wiggles" fit 
compared to the fit including the BAO feature. In addition, our cen- 
tral value and error on Dy/r^ are consistent when we fit ^0 only, or 
fo and ^2 simultaneously. 

Finally, we note that many of the small differences between 
the cosmological constraints presented here and those in our com- 
panion papers stem from slight differences in the best fitting value 



for Dv/r,(zdmE)/(-Dv/(Zeff)/r,(Zdrag))fiduciai- The Correlation function 
and power spectrum p ost-reconstruction "consensus" value from 
[Anderson et alJ ( 120121) is 1.033 ± 0.017; th is value was used in 
cosmological parameter studies in that paper. ISanchez et alj<2012h 
found 1.015±0.019, in agreement with the pre-reconstruction anal - 
ysis of the correlation function presented in And erson et alj ( l2012h . 
Once the WMAP7 prior on the underlying linear matter 
power spectrum p, is included and information from the Alcock- 
Paczynski effect is included through ^2, the standard ruler from 
the CMB allows us to infer (1 + Zcff)£'A(Zeff) and H{zcb) sepa- 
rately. Constraints from Model 3 (solid) and Model 4 (dotted) on 
(1 + Zcff)^A(Zcff) and /f(Zcff) are shown in Figures [8]and|9] Model 3 
further uses the WMAP7 ACDM prediction for /erg to disentangle 
the RSD and AP effects, and we find factors of 1 .3 and 1 .2 improve- 
ment in (1 -I- Zeff)£>/i and H errors when adopting this additional as- 
sumption: (1 + z,ff)DA = 2204 ± 44 (2190 ± 61) Mpc, H = 92.9+;^-5 
(92.4+;J-^) km s"' Mpc"' . In both models, the CMASS distance con- 
straints are consistent with what is inferred from WMAP7 alone in 
a ACDM cosmology: (l-FZeff)DA(Zeff) = 2ll3%H{z,ff) = 94.2!|;;]. 
We compare the two-dimensional constraints on (1 + z)Da and H 
from Models 3 and 4 with the prediction from WMAP7 for a flat 
ACDM model in Fig.[9l which shows that CMASS constraints on 
(1 + z)Da and H are only weakly correlated. 



6.4 Using our results 

Our results may be used to test cosmological models which share 
the assumptions we have adopted in this analysis. Most importantly, 
we have assumed adiabatic and scale-invariant primordial fluctua- 
tions, and that the transfer function was computed assuming the 
standard number of massless neutrino species, A^eff = 3.04. We have 
assumed that the linear growth is scale-independent, and account 
for non-linear corrections using perturbation theory within general 
relativity. The code to evaluate our theoretical prediction as a func- 
tion of the underlying linear matter power spectrum, cosmoxi2d 
is publicly availabl^j. For most purposes, however, our results can 



http://mwhite.berkeley.edu/CosmoXi2D 



© 0000 RAS, MNRAS 000,[T]tT] 



1 2 Reid et al. 



be well-approximated by the following multivariate Gaussian like- 
lihood for the parameters pj^, = [fcrf„F,{Dvlr,)l{Dvlr,)f,i^^i.^x\, 
which should be interpreted at Zeff = 0.57: 



Pid ■■ 



0.4298 
0.6771 
1.0227 



lO'C 



4.509868 2.435891 

2.435891 1.736087 

-0.01087251 -0.06287155 



-0.01087251 

-0.06287155 

0.3548373 



(17) 



(18) 



Deviations from the Gaussian likelihood are significant only for 
points > 3cr from the best fit values, where corrections to the like- 
lihood surface of ^0,2 will also become important. 



6.5 Key degeneracies and error budget 

In this section we examine the main sources of uncertainty in our 
measurements of fa f,(z^g) and -F(Zeff)- Both of these parameters af- 
fect the amplitude of the quadrupole, and so are partially degener- 
ate; Figure|9]shows that our measurements are sufficiently sensitive 
to distinguish the two through their differing scale-dependence. We 
will consider in turn the uncertainty due to the nuisance parame- 
ter cTp^Q describing the intrahalo velocities of satellite galaxies (the 
"finger-of-god effect"), uncertainty in the underlying linear matter 
power spectrum, the redshift-distance relation (in the case of pecu- 
liar velocities), and the peculiar velocity field (in the case of the AP 
effect). 

We first assess the impact of non-linearity in the covariance 
matrix on the eiTor budget. Taking the number of galaxies in the 
present analysis, assuming n = 3x lO"'* (h"' Mpc)~^, and using lin - 
ear theory to evaluate the Fisher matrix (as in lReid & Whitel201 ll) . 
we expect an uncertainty on fcrf, of « 0.021 when only bcr^ and 
faf, are freely varied. Using our mock covariance matrix, we find 
an uncertainty of 0.029 on /erg at fixed cTp^^ = 21 Mpc^ and for 
Fiin(fc), Dv, and F all fixed at their values in the mock catalog cos- 
mology. This ~ 40 per cent increase is primarily due to the non- 
linear/window function corrections to the covariance matrix high- 
lighted in Figure |5] If we instead use the mock based covariance 
matrix to fit for both Dy and F, fixing cTp^Q, P\m{k), and /erg, we 
find CT/r = 0.019. These results are summarized in Table|2] 



6.5.1 



Degeneracy with <Tp„Q 



At fixed P{k) and geometric parameters, the Fisher matrix analy- 
sis indicates a factor of 2 increase in the /erg error, to 0.042, when 
CTp^Q is marginalized over without any prior, compared to when it is 
fixed at CTp^Q = 21 Mpc^. The marginalized error on (J^^q in the for- 
mer case is 14 Mpc^. Therefore, the hard prior < cTp^Q < 40 Mpc" 
substantially reduces this source of uncertainty. Table|2]indicates an 
increase of only 15 per cent in the error on /erg, to 0.033, when we 
marginalize over cTp^^ within our hard prior. Similarly, marginaliz- 
ing over CTp^Q increases the error on F from 0.019 to 0.022. There- 
fore, further reduction in the uncertainty on (t\^q with more detailed 
modeling of the small-scale clustering would only allow a slight re- 
duction in the errors. However, since our fits indicate a slight pref- 
erence for cPpQQ = 40Mpc^ compared to the fiducial 21 Mpc^, in 
future work we will revisit our choice of prior after a re-analysis of 
small-scale CMASS clustering with a larger data set. 
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Figure 9. Upper panel: 68 and 95 per cent confidence regions for the co- 
moving angular diameter distance and expansion rate at z = 0.57 from 
CMASS anisotropic clustering constraints when /irg is varied over the 
WMAP7+GR flat ACDM prior (Model 3; solid) and when /erg is simul- 
taneously fit (Model 4; dotted). Lower panel: 68 and 95 per cent confidence 
regions for f<Tji(z = 0.57) and Alcock-Paczynski parameter F{z = 0.57) 
inferred from CMASS anisotropic clustering (Model 4; dotted). These two 
parameters are partially degenerate, and their differing scale-dependence al- 
lows us to constrain each separately. The solid contour shows the constraint 
when a WMAP7+GR flat ACDM prior is used on F. In both panels we 
show for comparison the predictions from WMAP7 when a standard GR, 
flat ACDM cosmology is assumed (dashed). 



6.5.2 Uncertainty in the underlying linear matter P(k) 

Of the three parameters defining the shape of the underlying lin- 
ear matter P(k) in our analysis, the uncertainty in oj^ = iich~ is 
the largest, particularly in the context of the constraining power of 
the CMASS measurements. For this comparison, we therefore hold 
Slijh- = 0.02258 and n, = 0.963 fixed, and examine how the central 
values and uncertainties in bag, fag or F change as a function of 
cjc, with all other parameters held fixed. Empirically, the uncertain- 
ties on these parameters do not depend much on ai^. Therefore, we 
can estimate the impact on the uncertainty in ai,. through the de- 
pendence of the central value of parameter p on 0)^., in units of the 
uncertainty on those quantities: 



Atoja^ 



(19) 
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Mai'ginalized parameters (Tjo-g ""fo-g <^F 



Model 



Fisher 


0.019 


0.021 


- 


- 


Fisher, <tI^^ 


0.024 


0.042 


- 


- 


- 


0.023 


0.029 


- 


- 


'^oG 


0.025 


0.033 


- 


- 


0.033 


0.033 


- 


- 


'^FoG' ''W' geometry 


0.033 


0.034 


- 


2 


- 


0.037 


- 


0.019 


- 


°"FoG 
4oG. PW 


0.038 


- 


0.022 


- 


0.046 


- 


0.022 


- 


o-g^Q, P(/:), growth 


0.046 


- 


0.026 


3 


'^oG 
<G. m 


0.047 


0.069 


0.042 


_ 


0.050 


0.069 


0.042 


4 



Table 2. We examine how uncertainty in various quantities entering our 
analysis impacts the 68 per cent confidence level intervals on parameters 
fecrg, fcTg, and F. For comparison, the first two rows show the predictions 
from a simple linear theory Fisher matrix analysis (as in Reid & White 
I2OI lb with « = 3 X 10""* (/i^'Mpc)"^^, the number of galaxies in the present 
analysis, and ignoring all window function effects. When <Tp^Q is marginal- 
ized over for the measurements (but not for the Fisher matrix analysis), 
we maintain a hard prior < o"|^q < 40 Mpc . Uncertainty in the under- 
lying linear P(k) is derived from WMAP7 data, under the assumptions of 
Gaussian, adiabatic, power law initial conditions with N^g = 3.04 mass- 
less neutrino species. Uncertainty in geometry [Dv(Zcff) and F(zcfi)] or 
growth [fcTg] is taken to be the uncertainty on these quantities derived from 
WMAP7 in a ACDM cosmology. The last column highlights the cases cor- 
responding to Models 2 though 4. 



where o",, is the uncertainty at fixed ate- For Gaussian probabil- 
ity distribution functions in P(p\a)c) and PiuJc), the uncertainty on 
p when marginalized over 0)^ is increased by Vl + s^. We find 
that 5 < 0. 1 for fag and s = 0.13 for F. Therefore, the current 
uncertainty in the underlying linear matter P{k) is negligible for 
the purpose of deriving constraints from the anisotropic cluster- 
ing of CMASS galaxies. This justifies the use of a fixed power 
spec trum shape in th e WiggleZ analyses of anisotropic cluster- 
ing JBlake. et al.ll201 Ih c). However, the uncertainty in P{k) does 
increase the uncertainty in bag, so applications suc h as galaxy- 
galaxy lensing/galaxy clustering combinations (e.g., iReves et al.l 
I2OIO.) should marginalize over this additional uncertainty. 



6.5.3 Uncertainty in ACDM geometry 

Table |2] indicates that marginalizing over the uncertainty in the 
geometric quantities (1 + ZciT)Z)A(Zeff) and i/(Zcff) [or equivalently, 
DviZeff) and F(zeff)] contributes negligibly to the uncertainty in fag 
for WMAP7 uncertainties when a ACDM redshift-distance relation 
is assumed. 



6.5.4 Uncertainty in ACDM growth of structure 

In the ACDM model, WMAP7 constrains Q,,, to ~ 1 1 per cent; this 
translates into a relatively large uncertainty in the predicted growth 
rate of structure (fag) of 5.5 per cent. Marginalizing over this un- 
certainty increases our uncertainty on F from 0.022 to 0.026. Weak 
lensing (e.g., Munshi et al. 2008) and/o r cluster abundances (e.g., 
IVikhlinin et Zll2009l : IRozo et al.ll2010h reduce the ACDM uncer- 
tainty on a combination of ag and Q.,„ similar to the OR prediction 
for peculiar velocities, crgnJJ;''', and thus could potentially be used 
to reduce the uncertainty on fag in the ACDM growth scenario. 



In the present work, we do not explore other dataset combinations 
besides WMAP7 and CMASS. 



6.6 Comparison with previous measurements atz~ 0.6 

The WiggleZ survey has recently analyzed the anisotropic cluster- 
ing of bright emissi on line galaxies over a broad redshift range 
JBlake. et al.l l2011bl la). Their growth rate constraints assumed 
a fixed underlying linear matter power spectrum and redshift- 
distance relation, for which our error is 0.033. With a factor ~ 4 
fewer galaxies, they achieve comparable precision and good agree- 
ment with our central value in their z = 0.6 bin: fag(z = 0.6) = 
0.43 ± 0.04. In their analysis they include modes up to fc,„ax = 
0.3/? Mpc"' and marginalize over a Lorentzian model to account 
for small scale nonlinearities. We find that the quoted WiggleZ er- 
ror is in good agreement with our Fisher matrix prediction if we as- 
sume n = 2.4x lO"'* (/^"'Mpc)"^ b = l.l,Ng,,]= 60227, a Gaussian 
damping croauss = 300 kms"' and 5n,i„ = l.l5Kk^\^ - 12/j"'Mpc, 
i.e., a factor of two smaller scale than we have adopted for our 
analysis. Exploration using the Fisher matrix suggests that the dif- 
ference in number densities between CMASS and WiggleZ has 
a negligible impact on the uncertainties, while the lower bias of 
their sample implies a 15 per cent (10 per cent) improvement at 
•^'rain = 25 (12)/;"' Mpc. By far the dominant difference arises be- 
cause they are fitting out to k^^ = 0.3 h Mpc"' ; their error would 
increase by a factor of 2 if they adopted our minimum scale. An- 
other issue is that our mock covariance matrix, which accounts for 
both nonlinear growth of structure and our complex survey geome- 
try, yield errors 40 per cent larger than our naive Fisher matrix anal- 
ysis would predict, with the difference between the linear and mock 
covariance matrices increasing on small scales (for ^2)- While it is 
not clear how this difference scales with galaxy bias or depends on 
survey geometry, the WiggleZ use of linear theory covariance ma- 
trices down to much smaller scales could cause their uncertainty to 
be underestimated. 

Another recent WiggleZ analysis ( IB lake, et al.ll201 id) simul- 
taneously fits F = 0.68 ± 0.06 and fag = 0.37 ± 0.08 at z = 0.6 
for a fixed underlying P{k), and makes use of modes with k < 
0.2/? Mpc"', which can be compared with the last two lines of Ta- 
ble [2] and Model 4 in Table [T] As in the previous discussion, our 
best fit values for these parameters are consistent, and using clus- 
tering information on small scales (i.e., larger k^^^) permits a factor 
of ~ 1 .4 tighter constraint on F at a fixed number of galaxy spectra 
compared with our clustering analysis of CMASS galaxies. How- 
ever, the higher bias of CMASS galaxies permits relatively tighter 
constraints on F than fag as compared with WiggleZ. 

The WiggleZ and BOSS CMASS surveys target very different 
galaxy types, which will have different nonlinear properties and 
modeling uncertainties. Without a detailed study it is impossible to 
understand the robustness of various assumptions for the nonlin- 
ear distortions, particularly when multiple sources of nonlinearity 
contribute. These can potentially cancel when examining cluster- 
ing, meaning that goodness-of-fit must be used carefully when as- 
sessing robustness and potential biases in a model fit. In this work 
we have adopted a conservative approach based on the best model 
available that has been compared with a large volume of W-body 
mock galaxy catalogs, and only fit over those scales where we are 
confident that the signal is dominated by the quasilinear velocity 
field of interest, and where the impact of small-scale random mo- 
tions can be simply modeled and marginalized over. 
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Figure 10. Joint constraints on n^ from tlie combination of CMB distance 
and normalization priors and our measurements of Da , H, and /erg . Tliese 
constraints are independent of the growtli and expansion liistoiy at redsliifts 
lower than our galaxy sample, but assume a flat ACDM cosmology between 
the CMASS sample and the CMB, as well as the other model assumptions 
detailed in Section ISTI 



7 DISCUSSION 

We have analyzed the anisotropic clustering of BOSS DR9 CMASS 
galaxies with an accurate analytic model for the monopole and 
quadrupole correlation functions, which we have validated using 
a large volume of A'-body based mock galaxy catalogs. The com- 
bination of the BAO standard ruler and the Alcock-Paczynski ef- 
fect allows us to separately constrain the comoving angular diam- 
eter distance and Hubble expansion rate at the effective redshift 
of our sample, z = 0.57, while redshift space distortions allow 
us to constrain the amplitude of the peculiar velocity field, a di- 
rect measurement of the growth rate of structure, dcrg/dlna. Ta- 
ble[T]summarizes our constraints under several assumptions. Figure 
|9] shows the degeneracy between the Alcock-Paczynski parameter 
F and the growth rate of structure in our measurements; this ex- 
plains why our constraints improve considerably if we make further 
model assumptions about the geometry or growth. In the most gen- 
eral case where all three parameters vary independently, we find 
do-s/dlna = 0.43!!!i!g, (1 + Zeff)£»A(Zetf) = 2190 ± 61 Mpc and 



H{zcb) = 92.4!:! 5 km s"' Mpc"'. 

To illustrate the cosmological constraining power of our 
measurements, we summarize our results as three distinct tests 
of a minimal ACDM cosmology (see Section [STt that link the 
observed CMB anisotropics at z a 1091 and CMASS galaxy 
fluctuations at z a 0.57, independent of the expansion history 
and growth of structure in the universe at redshifts below our 
sample (z ;£ 0.57). We use a flat ACDM model to illustrate our 
constraining power on the behavior of the universe at z > 0.57, 
and e xplore more general co smological models in our companion 
paper, ISamushia et al.l ( l2012h . 

Is the shape of the power spectrum of matter density fluctua- 
tions inferred from the CMB consistent with the one inferred from 
galaxy fluctuations after a factor of ~ 4 X 10' ampliflcation? 



ACDM evolution and galaxy biasing (~ 10 per cent on the scales 
we analyse), the shapes of the linear matter power spectrum 
inferred from the CMB and galaxy clustering are consistent: our 
best fit ACDM model gives ;jf^ = 39.3 for 44 degrees of freedom. 
We quantify this statement further using our fit to the location 
of the broad turnover in P(k), which indicates the horizon size 
at matter-radiation equality and thus the physical matter density, 
assuming the radiation density is known. Combining our constraint 
(Q.,„h- = 0.142 ± 0.01 1) with those in the literature for the SDSS-II 
LRG sample and WiggleZ yields Q.,X = 0.137 ± 0.0064, where 
the error neglects the expected small but non-zero covariance 
between the galaxy samples, and also fixes Q^/i" and n, to best fit 
CMB values. The CMB constraint is only slightly more precise 
than the combined galaxy measurement {Q.„,h^ = 0.1 334 ±0.0056), 
and the two are in excellent agreement. With this test passed, 
in the rest of our analysis we make use of the full matter power 
spectrum (rather than just the BAO feature) as a standard ruler 
in galaxy clustering measurements. See ISanchez et al.l ( 120120 for 
an exploration of other cosmological models using the CMASS 
monopole correlation function shape. In particular, the good agree- 
ment between constraints on the linear matter power spectrum 
from the CMB and CMASS galaxy clustering limit the allowed 
contribution from species such as massive neutrinos which induce 
scale-dependent growth. 

Do our constraints on the geometry of the universe require 
dark energy at z<i Q.51? 

Distance constraints from both the CMB and the BAO fea- 
ture are determined relative to the sound horizon at z* « 1091, 
the redshift of decoupling, and at Zdrag ~ 1020, the redshift when 
baryons were released from the Compton drag of the photon^j. 
Within the constraints of our cosmological model assumptions (see 
Section [STt . r^izt) and rj(zdrag) depend only on Q.,„h^ and Q.),h^. 
We project the observed difference in comoving angular diameter 
distance AZ)^ = (1 + 1091)Da(z = 1091) - (1 + 0.57)Oa(z = 0.57), 



A£»4 



H{z) 



(20) 



onto the parameter Q.\ in a flat ACDM cosmology and marginalize 
over the other parameters. The result is shown by the dashed 
line in Figure [TO] Qa = 0.76 ± 0.04. We quote the maximum 
likelihood and 68 per cent confidence region around it unless 
otherwise noted. We have purposely chosen a variable that is 
independent of the expansion history at z < 0.57. Similarly, the 
CMB constraint on il^h^ with our measurement of Hiz = 0.57) 



requires Da = 0.71^ 
curve in FigurefTOt. 



between 



0.57 and z = 1091 (dotted 



Does the observed growth rate of structure at z = 0.57 require 
dark energy ? 

Because WMAP places such a tight constraint on the ampli- 
tude of curvature perturbations deep in the matter-dominated epoch 
at fc = 0.027 Mpc"' (1.8 per cent), we can translate our measure- 
ment of /(Tg into a constraint on dZ)/d In a, where D is the usual 
growth function that depends only on fi,„ or Qa in a flat ACDM 



If the dominant component of the energy density is "cold", 
then in the linear regime perturbation growth in General Relativity 
is scale-independent. Modulo our corrections for non-linear 



^ We follow the WMAP7 analysis and adopt the fitting formulae for z* 
and zji-ag from .Hu & Su gi yama , (.1996.) and. Eisenstein & Hu, y_99^, respec- 
tively. 
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cosmological model: 



D(z)^ 



5n„ H(z) 



■■m r 

Ho J 



I/(l+z) 



da Hi 
[a//(a)]3 



(21) 



To do so, we require the additional cosmological model assump- 
tions listed in Section 15.11 We marginalize over the WMAP7 
uncertainties in the parameters that convert curvature perturbations 
to the integrated amplitude of matter perturbations on scales of 
R = 8/0.7 = 11.4 Mpc, (Tr, namely Q.„,h-, Q.i,h^, and «,. Section 
I5.3l details the relation between gtr and erg, which can be regarded 
as equal for our purposes. This comparison between the fluctuation 
amplitude at z = 1091 and z = 0.57 requires Qa within [0.59, 0.81] 
(central 68 per cent confidence), with a maximum likelihood at 
0.76. The distribution is shown as the dot-dashed curve in Figure 

m 

In combination with the CMB, both our geometric and growth 
rate constraints require a value of Q,\ at 0.57 < z < 1091 that is 
consistent with the concordance model expectation. Combining 
all three constraints with the WMAP7 distance and normal- 
ization priors, and marginalizing over Q.bh- and Q.,„h^, we find 
n^ = 0.74°-^j,'|5 (solid curve in Figure [TO). 

Our analysis requir ed further assumption s compared with 
a BAO-only analysis l iAnderson et al.l 120120 that constrains 
Dv{Zcft)l''s{Zirig)- However, we were able to perform tests of the 
scale-independence and growth rate of cosmic structure between 
recombination and z ~ 0.57, as well as to break the degener- 
acy between (1 -t- Zcfd^A and H using the Alcock-Paczynski test. 
So far, our measurements do not unveil any deviations from the 
minimal ACDM model we have examined. Under the assump- 
tion of a ACDM cosmology extending to z = 0, the BAO feature 
adds the most constraining power to WMAP7 on ACDM param- 
eters; we find n„,/i^ = 0.136 ± 0.0035, n„, = 0.283 ± 0.017, and 
Ho = 69.3 ± 1.5 km s"' Mpc"' in this model. We anticipate statisti- 
cal improvements on these results with the completed BOSS galaxy 
dataset covering a footprint three times larger, as well as develop- 
ments in the theoretical modeling that will allow tighter cosmolog- 
ical constraints. 
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APPENDIX A: OBSERVATIONAL UNCERTAINTIES 
A I Uncertainties in radial distribution 

To compute overdensities of the galaxy field, defined as 



S{z,a) 



pq(z, n) - p(z, fl) 
po(z, n) 



(Al) 



we must know the unperturbed density field po(z, f2). While the an- 
gular selection function of the survey is usually well known, the 
radial distribution is not easy to model accurately. Usually, the un- 
perturbed radial distribution of galaxies is modeled from the dis- 
tribution of observed redshifts by either shuffling them or splining 
with a smooth curve. Different ways of constructing a random cat- 
alog W2ll_^esultjn_different_£stin^^ (see, 
e.g.. lSamushia. Percival, & Raccaneiiill20Ill : lRoss et al.ll2012h . 

To estimate the magnitude of this effect on the measurements 
of the moments of correlation function we take 600 mock catalogs 
of CMASS sample and apply simplified version of our analysis to 
the measurements produced using different ways of reconstructing 
n(z). We reconstruct n(z) first by shuffling "observed" redshifts in 
mock catalogs and then by splining that distribution with 10, 20 
and 30 node cubic spline fits. 

We find that the ways of reconstructing n(z) for random cata- 
logs do not affect the RSD measurements significantly. In our cur- 
rent analysis we are using the "shuffled" catalogs since they intro- 
duce the least bias in the measurements of ^/(r) (.Ross et al.2OI20 . 



A2 Effects of Close Pairs, redshift errors and redshift 
failures 

We are unable to obtain redshifts for ~ 5 per cent of the galaxies 
due to fiber collisions - no two fibers on any given observation can 
be placed closer than 62". At z ^ 0.5 this 62" exclusion corre- 
sponds to 0.4 /i"' Mpc. Redshifts for some of the collided galaxies 
can be reclaimed in regions where plates overlap, but the remain- 
ing exclusion must be accounted for. We account for fiber-collided 
galaxies by assigning its weight to it s nearest neighbor on the sky . 
Tests on mock catalogs presented in JGuo, Zehavi. & Zhendl201lh 
indicate that the nearest neighbor correction adopted in this work 
is accurate to better than 1 per cent for both ^o and ^2 at the scales 

used in our analysis. 

Redshift failures are discussed in detail in lRoss et al.l ( l20I2h : 
we also correct for them with a nearest-neighbor upweighting 
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scheme. Redshift measurement errors smooth the apparent galaxy 
density field, in the same fashion as described with our nuisance 
parameter cTp^^. The median redshift error for our sample is 42 
km s"', which translates into an additive contribution to crp„e < 1 
Mpc^. 



APPENDIX B: ACCURACY AND IMPLEMENTATION OF 
THE THEORETICAL MODEL 

Bl Wide Angle Effects 

Eauationsll0land[T2l assume that a "plane-parallel" approximation, 
which states that the sample is far enough from the observer so 
that all line-of- sights are parallel to each other, is accurate enough. 
Thi s approximation will f ail at some scale for wide surveys (see 
e.g. IPapai & Szapudill200 8). We approximate the scale dependent 
m agnitude of wide-angle effects for our sam ple in a similar manner 
to lSamushia. Percival, & Raccanellil ( 120111) . The effective redshift 
of our galaxy sample is z = 0.57, which in the best-fit WMAP7 
cosmology corresponds to a comoving distance of approximately 
1500 /j^'Mpc. The largest scale that we consider in this analysis 
is 160 /i^'Mpc, which corresponds to an opening angle of about 3 
degrees. We estimate the wide-angle corrections to be at most 10 
per cent of the statistical eiTors on the largest scales and do not try 
to correct for this effect in current analysis. 



B2 cosmoxild code Implementation 

iReid & White! j201lh demonstrated that Eqn. [12] provides an 
accurate description of ^iira, ''ti) when real space clustering and 
velocity statistics inside the integrand are measured directly from 
A'-body simulations. Moreover, in the regime where b k: 2 (appro- 
priate for CMASS galaxies), the real space clustering and velocity 
statistics can be computed analytically with sufficient precision to 
predict ^lir^, r^). In detail, this "sweet spot" in the precision of the 
model arises because the functions V\2{r) and a^^^^.J^r) entering 
Eqn. [T2] were evaluated as a function of cosmological parameters 
in standard perturbation theory accounting for only the linear bias 
of the tracer. At the redshift of interest, the second order bias bit 
crosses zero near 1 -I- b\L = 2, and the calculation of these functions 
neglecting b2L is sufficient for our purposes. We therefore caution 
against the use of our code for tracers with bias s ubstantially 
differe nt from 2. Note that ^'g(r) is evaluated in LPT JMatsubaral 
uOOSbf ) and includes the contribution to the real space clustering 
from second-order Lagrangian bia^ 



Standard perturbation theory is known to have inaccuracies 
in describing the BAG feature, and we found that on BAG scales 
errors in our perturbation theory calculation of pairwise velocity 
statistics caused inaccuracies in the prediction for ^lir^, r„). How- 
ever, on scales s > 70 /;"' Mpc, the redshift space version of LPT 
I MatsubaraluOOSbh is very accurate (see figure 2 of iReid & White! 
!201l!) . We therefore interpolate between our evaluation of Ean.ll2l 
at smaller scales and LPT on large scales, with the transition fixed 
at 100 Mpc. 

Finally, we model the effect of galaxy intrahalo velocities (tra- 
ditional "fingers-of-god") by convolving our model f (^p, r„) with 



^ We relate the second order Lagrangian bias, b2L, to the first order bias, 
b[[_, through the peak background split; these parameters are not varied in- 



dependently. 



an additional Gaussian velocity dispersion of variance (Tp^Q. Using 
the Gaussian form allows o"p„Q to be included directly in the Gaus- 
sian in Eqn. [T2|for faster evaluation of our model as we explore 
cosmological parameter space. 

We have developed a code to numerically evaluate Eqn. [T2| 
as well as all the relevant perturbation theory integrals as a func- 
tion of an input linear matter power spectrum and nuisance param- 
eters b and (Tp^Q. The Alcock-Paczynski effect is easily incorpo- 
rated using Eqn. [13] before computing Legendre polynomial mo- 
ments from ^lira-, r„}. The internal units of cosmoxild are Mpc, in 
which the underlying linear matter power spectrum is most tightly 
constrained. The code is publicly availablqS 



B3 Model Accuracy 

Because we have such a large volume of simulations, we can use 
the difference between the theoretical model at the known cosmo- 
logical parameters of the W-body simulation and the measured cor- 
relation function from the mock galaxy catalogs to quantify our 
theoretical systematic error; we find Ax^ = 0.29 at the best fit value 
of cTp^jj. We also compare the scale dependence of the model error 
with the 5 parameters we are fitting, and find < 0.25cr shifts com- 
pared to the unmarginalized uncertainties on all parameters (i.e., 
the uncertainty on each parameter if all the others were known per- 
fectly). Therefore we conclude that our systematic error is negligi- 
ble in the context of this analysis. 

Another concern is that the model becomes inaccurate rapidly 
on scales smaller than our minimum fitting scale, .s,ni„. When 
Alcock-Paczynski parameters vary, scales smaller than iniin con- 
tribute to the model. However, we have verified that for variations 
of 10 per cent in Dy and/or 20 per cent in F (i.e., much larger than 
the final uncertainties), the theoretical error induces < 0.5cr shifts 
in all parameters compared with the unmarginalized uncertainties. 

Finally, we point out that unlike in linear theory, our model 
depends on b, f, and erg separately, rather than the only the combi- 
nations bcTg and /erg. However, at the particular berg and fcrg val- 
ues of our sample, we find that changes in the predicted ^0.2 with 
(Tg can be absorbed by changes in the value of berg. Quantitatively, 
for a ±10 per cent change in away from our fiducial cr8(ZefF) = 0.61 
and at fixed P{k), AP parameters, and cr^^Q = 21 Mpc^, bag shifts 
by a; ±1.3cr with no measurable shift in the central value of fag. 



B4 Propagating uncertainties in the galaxy-halo mapping to 

0"foG 

As shown in Figure |6] our model with o"p^Q = 21 Mpc" fits the 
(o ijs) of our mock g alaxy catalogs, based on the best-fitting HOD 
in !White et al.! ( !201 ID . In this section we quantify how uncertainties 
in both the theoretical modeling and data analysis cause uncertainty 
in the expected value of o"p„Q for the CMASS sample. We address 
several aspects of this problem separately. 

• One-halo vs two-halo contributions to fo,2 and fiber collision 
corrections: The formalism of the halo model distinguishes be- 
tween "one-halo" and "two-halo" pairs depending on whether the 
two galaxies occupy the same or different halos. In Figure IbTI we 
show the total change in fo,2 due to satellite galaxy intrahalo ve- 
locities (IHV) by the dashed (dotted) curves in our mock catalogs. 
The dash-dot and solid curves show the contributions to this change 



"* http://mwhite.berkeley.edu/CosmoXi2D 
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Figure Bl. The dashed (dotted) hnes show the fractional change in fo 
(^2) in our mock galaxy catalogs due to satellite galaxy intrahalo veloci- 
ties (IHV; also shown in Figure |6). We isolate the contribution of pairs of 
galaxies occupying the same halo in our mock catalogs (i.e., "fingers-of- 
god"), shown as the dot-dashed (solid) curves for fo (ft)- On the scales of 
interest, the dominant effect of IHV is a net diffusion of pairs from small 
scales (where f is larger) to larger scales. 



in ^0,2 from one-halo pairs, which are locahzed to along the LOS 
with ru < I h^'Mpc. Because the one-halo pairs contribute to such 
a small rf/j, they can be neglected on the scales used in our cos- 
mological parameter fits, s > 25 A^'Mpc. This fact is important 
to establish since our method for fiber collision corrections will 
correctly recover the distribution of pair separations for pairs of 
galaxies with separations larger than the fiber collision scale, but 
suppresses the contribution of pairs of galaxies at the fiber colli- 
sion scale (ro- < 0.5 /i^'Mpc). 

• Uncertainty in the HOD at fixed cosmology: Uncertainties in 
the HOD parameters will introduce an uncertainty in cr^^^- While 
CTp^Q is roughly proportional to the satellite fraction, it also depends 
on the distribution of host halo mass - increasing a and k increases 
host halo mass at fixed s atellite fraction, wh ich increases o"p„Q. We 
use MCMC chains from I White et alj ( 120 111) to estimate our uncer- 
tainty on o"p„Q at fixed cosmology to be 6 Mpc^. 

• Breakdown of the "central" galaxy assumption: The analysis 
of ISkibba et alj ( 120111) suggests that the brightest galaxy in a halo 
is not always the "central" one. We test the impact of relaxing our 
assumption that the velocity of mock central galaxies have no intra- 
halo velocities by assigning them the intrahalo velocity of a random 
dark matter particle halo member in our simulations some fraction 
p of the time. When we have more than one galaxy in a halo, we 
assume that the chance of not including the "central" galaxy in our 
sample is lower, oc p"s"'. For p = 0.3, we find a-p^^ at our fiducial 
HOD is increased by 9 Mpc^ . 

• Variations in the halo mass function with cosmological pa- 
rameters: A broad range of observations shows good agreement 
between the concordance ACDM halo mass fun ction and the mul- 
tiplicity of galaxy groups an d clusters (e.g., IRozo et al.ll20ld : 
lAllen. Evrard. & Mantzll201 11) , so this uncertainty is subdominant: 
even if our fiducial HOD masses were scaled by a factor of 2, 
cr^^Q oc M"'^ would change by 1 1 Mpc" . 

Given the above considerations, we adopted a generous hard prior 



APPENDIX C: MARKOV CHAIN MONTE CARLO 
(MCMC) METHODS 

We adopt a hybrid MCMC/importance resampling approach to ex- 
plore the BOSS likelihood surface with various priors imposed 
from the WMAP7 likelihood in Models 1-4. This approach is nec- 
essary in our case because our model evaluation is slow, and we 
must marginalize over /)0"g and cTp^Q at each point in parameter 
space that we consider. This section describes our methods in each 
case. 



CI Importance resampling 

WMAP7 MCMC chains are publicly availablfl These chains 
provide a fair sample of the WMAP7 likelihood s urface . 
Importance resamplin g dGilks, Richardson & Spiegelhalteij|l99q : 
iLewis & Bridlell2002l) allows us to compute how constraints on the 
model parameters change given an additional constraint by eval- 
uating the new likelihood at a subsample of the original MCMC 
chain, multiplying the original weight of each element by the new 
likelihood, and then recomputing confidence intervals. 



C2 Model 1: WMAP7+CMASS ACDM 

For this model we marginalize only over two parameters, so it is 
feasible to directly compute a marginalized CMASS likelihood: 



(^A 



CMASSl.t'ACDM 



,) 



/ 



dbo-g do-^^Q e 



-J^CMASS<''ACDM.*'rg,tr2^,G)/2^ ^Q^^ 



on cTp^Q between and 40 Mpc 



We then use importance resampling of the WMAP7 chain parame- 
ters. 



C3 Model 2 WMAP7-I-CMASS ACDM geometry, free 

growth; Model 3 WMAP7-I-CMASS ACDM growth, free 
geometry 

In Models 2 and 3, for each linear matter power spectrum, we 
must vary three or four extra parameters describing the galaxy 
clustering observations. We therefore explore the likelihood sur- 
face at each point in the WMAP7 chain by MCMC, and thereby 
sample the CMASS likelihood distribution -PC^cmassI^acdm)- In 
Model 2, ScMASs = {ftcrg,cr^^|-.,/crg), and in Model 3, ^cmass = 
{fccTg, cTp^Q, Dv, F}. MCMC chains at a fixed (?acdm sample param- 
eter space proportional to /'(0cmass|6'acdm)> but in order to compute 
the marginalized likelihood of fcTf, in Model 2 or Dy, F in Model 
3, we must determine the relative likelihood of the MCMC chains 
evaluated at different ^acdm- Since 

yVg-^^(ecMASS.«ACDM)/2 = /'(^cMASS , ^acdm) = P{ecM^&^\eKCDM)P{0KCl>M) 

(C2) 

where A' is an overall normalization, we can combine the x^ 
computed at any point with the probability density estimated by 
our MCMC to determine the relative normalization of -P(0acdm)- 
In practice, we first normalize each MCMC distribution so that 
/^(ScmassI^acdm) = 1, find e^^^^^iO^cDu) with the minimum 
value of ^^ in each chain, and integrate P(6cmass|Sacdm) in a small, 
fixed size region of parameter space around Gquass' which we 



http:/Aambda. gsfc.nasa.gov/product/map/dr4/parameters.cfm 
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call /)(0^^^gg,rf6cMASs)-L|The relative weight of each point in the 
WMAP-only MCMC chain is then determined by the CMASS-only 
likelihood value and p: 



wC^acdm) = e 



-x-(e* 



(P(^CMASS''^^C 



s)r 



(C3) 



We find that our constraints are the same if we neglect this volume 
weighting factor p"' , indicating that the effective volume of param- 
eter space allowed by the CMASS measurements does not strongly 
depend on the underlying cosmological parameters 0acdm when ex- 
ploring the region of this parameter space allowed by WMAP7. 

C4 Model 4 WMAP7-I-CMASS, free growth, free geometry 

In this case, we use WMAP data to provide a prior on the shape of 
the linear matter power spectrum, which is well approximated by 
a multivariate Gaussian in the parameters Slch', Clblr^, iig. WMAP7 
constraints on these parameters primarily come from ratios of peak 
heights and the overall shape, rather than the locations of the peaks 
(which are sensitive to the angular diameter distance to the last scat- 
tering surface). Therefore, marginalized likelihood for these param- 
eters is nearly independent of the adopted model for the low red- 
shift expansion history (i.e., whether Q.^, or w are freely varied). 
Thus we are able to make use of the CMB information on the un- 
derlying linear matter power spectrum that is independent of the 
model for the low redshift expansion history, and does not contain 
information on the distance to the last scattering surface. By us- 
ing the linear matter power spectrum as a "standard ruler", we are 
able to infer information about the geometric parameters Dv(Zcti) 
and F. In this case we run an MCMC chain with the following 8 
parameters, adding an additional multivariate Gaussian likelihood 
representing the CMB prior on the three linear matter power spec- 
trum parameters: {nch^,Q.iJi^,n„ Dv(Zsff), F, b<rg,fcrg,cr^^]. 
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